Brush-like Interactions between Thermoresponsive Microgel Particles 
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Using a simplified microstructural picture we show that interactions between thermosensitive 
microgel particles can be described by a polymer brushlike corona decorating the dense core. The 
softness of the potential is set by the relative thickness Lo of the compliant corona with respect 
to the overall size of the swollen particle R. The elastic modulus in quenched solid phases derived 
from the potential is found to be in excellent agreement with diffusing wave spectroscopy data and 
mechanical rheometry. Our model thus provides design rules for the microgel architecture and opens 
a route to tailor rheological properties of pasty materials. 
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Thermoresponsive microgel particles are a hybrid be- 
tween a colloidal and a polymeric system with proper- 
ties that can be tuned externally [2H4]. Most of the 
previously studied stimuli responsive microgel systems 
are based on poly(N-isopropyl-acrylamide) (PNIPAM), a 
polymer which has a lower critical solution temperature 
(LCST) of approximately 33°C 0. Above the LCST, 
the microgel particles expel water and are collapsed. The 
typical size of a collapsed microgel particle is in the range 
of 0.2-1 (tm. Upon lowering the temperature the parti- 
cles swell to about twice their original size. Responsive 
microgels thus provide the possibility to fabricate smart 
colloidal materials for applications as viscosity modifiers, 
carrier systems, bioseparators, optical switches or sensors 
d, 0, Moreover, due to their tunability they are ideal 
model systems to study the phase behavior, glass transi- 
tion and jamming in dense colloidal dispersions @, 0]- 

Particles come into contact and form a viscoelastic 
paste below the LCST if the polymer density inside the 
swollen particles approaches the total polymer density. A 
number of rheological studies have revealed the apparent 
divergence of the viscosity at the transition point and the 
emergence of an elastic shear modulus [3j, |j| . Since the 
particles consist of polymers and solvent in equilibrium 
they are elastically compliant, and thus elastic moduli do 
not diverge at the jamming transition, unlike the behav- 
ior of rigid colloidal particles [1, 0] . Several experimen- 
tal studies show that in this regime the bulk modulus 
as a function of the effective volume fraction scales as a 
power law G p oc ip 1 ^?^ with values of n ranging from 
n = 9 to n = 22 [3|, These findings indicate that 
the particle interaction potential ip(r) oc r~ n also fol- 
lows a power law, at least over some significant range of 
length scales [H, 0] • Unfortunately a detailed model for 
the origin of these interactions has been missing. This 
is mainly due to the fact that modeling the interaction 
between swollen particles is complicated by the heteroge- 



neous microstructure arising from the faster reaction rate 
of cross- linking compared to polymerization flpj j . In this 
Letter we show that interactions between thermosenstive 
microgel particles can be modeled by apolymer brushlike 
corona decorating the dense core [Ill4l3j . The softness of 
the potential is set by the relative thickness of the compli- 
ant corona with respect to the overall size of the swollen 
particle. Our PNIPAM microgel particles have been syn- 
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FIG. 1: PNIPAM hydrodynamic radius Rh as obtained from 
dynamic light scattering. Dash-dotted line shows the inter- 
polation curve used to calculated the effective volume frac- 
tion $ e // oc Rjj. The sketch (drawn to scale) illustrates the 
swelling of the particles: Our experiments suggest that the 
particles consist of a highly cross-linked core with radius Rb 
and a corona of uncross-linked polymer chains that can swell 
to a brush of thickness Lo — R — Rb under good solvent con- 
ditions (left). 

thesized by standard methods and characterized as de- 
scribed in Ref. 14]. The high cross- linker to monomer 



ratio of 5.3 mol % N,N'-methylene-bis-acrylamide (BIS) 
leads to relatively rigid microgel cores with, on average, 
one molecule of cross-linker per 19 molecules of NIPAM 
monomer. In Fig. [T] we show the temperature depen- 
dence of the hydrodynamic radius Rh determined by 
dynamic light scattering. The particle radius in the col- 
lapsed state is Ro — 250nm for which the internal poly- 
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mer volume fraction is approximately (j>o w 70% 14 1. 
For T= 20° C the particles reach a hydrodynamic radius 
of Rh — 440nm, almost doubling their size. As parti- 
cles swell they acquire a radially inhomogeneous density 
profile. While the exact density distribution <f>(p) in the 
swollen state remains unknown, scattering experiments 
clearly indicate the presence of a moderately swollen core 
and a highly swollen corona 0, ll4T - [l6j ] . For our particles 
we find from static light scattering a polymer density of 
4> ~ 35% in the core and <fr < 10% in the corona (at 
p > 0.8R H ) 0. Classical Flory-Rehner theory 0, [}J§ 
for swollen polymer gels provides estimates for the degree 
of polymerization between cross-links N c in the core and 
corona regions. Using literature data for the statistical 
segment length a ~ 0.81nm [l8| and the Flory interac- 
tion parameter x El of NIPAM, we find N c ~ 10 in the 
core of the particle, in agreement with the stoichiomctry 
of the microgel synthesis. Moreover, for typical corona 
densities we find N c 3> 100, suggesting an essentially 
uncrosslinked corona region for our polymerization con- 
ditions. Below the LCST ~ these corona chains extend 
into the good solvent and form a swollen brushlike layer 
of thickness Lq. While the transition between the core 
and corona regions is not an abrupt one, we adopt the 
simplified scaling picture of Alexander ll| of the brush 
as a collection of identical chains of length N grafted at 
moderate-to- high surface density to a reference surface at 
p = Rb- In good solvent conditions, the thickness of such 
a brush scales as L ~ As~ 2 / 3 a 5 / 3 , where s ~ aa^ 3 ^ 4 is 
the average separation between grafting sites. Note that 
other more realistic models also predict Lq ~ N [l2| . so 
the simplicity of our model should not adversely affect the 
results at the scaling level. Results from static light scat- 
tering indicate typical values of 4>b ~ 3% and Lq ~ lOOnm 
for our particles [14J . From these values we can estimate 
s ~ lOnm and N ~ 700. This means the corona density 
is well above the overlap concentration <j)* ~ N~ 4 / 5 and 
s < Lq which shows the consistency of our arguments. 
We note that due to the high degree of cross-linking the 
stiffness of the core region is substantially higher than 
the soft brushlike corona Uj . More generally, as long as 
the number of monomers between cross-links is less than 
the number of monomers in a corona chain, the core will 
be less compliant than the corona region. As we are 
well within this regime, we may treat the core as an ef- 
fectively incompressible solid for the weak to moderate 
particle deformations considered. Here, the choice of a 
reference surface, Rb, is somewhat arbitrary but should 
reflect the boundary between core and corona behavior. 

We now turn our attention to the properties of a dense 
microgel assembly with a total polymer concentration of 
about 135 mg/ml [14]. From laser scanning confocal 
microscopy studies of a diluted sample we determined 
the particle number density of the dense system to be 
r\ = 2.24/yum 3 ± 2%. Over the range of temperatures 
covered, the effective particle volume fraction 3> e // thus 



increases from about 0.2 to 0.9. We use two-cell diffus- 
ing wave spectroscopy (DWS) in the transmission geom- 
etry to study particle diffusion, dynamical arrest, and 
build up of elasticity [20(. DWS, dynamic light scatter- 
ing in the highly multiple-scattering limit, provides in- 
formation on the thermally-induced fluctuations of the 
microgel particles on the nanoscale via the temporal in- 
tensity autocorrelation function g^ (t) of multiply scat- 
tered light. In particular, the ensemble-averaged mean 
square displacement (MSD), (Ar 2 (i)), can be extracted 
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from gi (t) . Details of the setup can be found in Ref. 
The sample is kept in a glass cuvette with inner dimen- 
sions 10 x 2 mm (Hcllma, Germany). To avoid crystal- 
lization the sample is loaded hot, quenched below 20°C, 
and subsequently the temperature is increased in steps 
of 0.5°C or 1°C. At each temperature the sample is kept 
for about lOmin to equilibrate and ga(t) is recorded over 
an additional 10 min. For our microgel suspensions, g2(t) 
(partially) relaxes in a characteristic time r that depends 
on temperature. In the case of more dilute suspensions, 
g2(t) exhibits a roughly exponential decay characteristic 
of simple diffusion, with (Ar 2 (i)) = 6D s (4>)t. However, 
for concentrated suspensions, the behavior is more com- 
plex. At temperatures between 27°C and 30°C, g%{t) 
exhibits a slightly stretched exponential decay and a cor- 
responding power-law MSD, (Ar 2 (t)) ~ t p with /3 <, 1, 
characteristic of a crowded suspension of soft, repulsive 
particles. 




FIG. 2: Thermally activated nanoscale motion. Mean square 
displacement (Ar 2 (£)) of the microgel particles from diffusing 
wave spectroscopy. 

However, at T = 27°C and below the particles have 
swollen into contact and the short-time decay of gi (t) is 
only a partial one, leading to a plateau value of the MSD 
characteristic of the elastic properties of the jammed sus- 
pension FigureCHshows plots of (Ar 2 (i)) versus time 
t for a range of temperatures above and below the jam- 
ming transition. 

We first analyze the influence of crowding on the tem- 
perature dependence of the self diffusion coefficient D s ((f)) 
for T > 30°C. As shown in Fig. [3] our results are 
well described by the semiempirical expression of Lion- 
berger and Russel for hard spheres 23[: Do/D s ($> e ff) = 
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[(1 - 1.56 $ e //)(l -0.27 $ e// )] \ where $ e// = rj x 
47ri? 3 /3. Here we have introduced an effective steric ra- 
dius R = Rh x 1.027 that reflects quantitatively the hard- 
sphere like behavior in the liquid state. At tempera- 
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FIG. 3: Microgel liquid-solid transition as a function of tem- 
perature. The transition is driven by the change in effective 
particle volume due to swelling. The reciprocal short-time 
diffusion coefficient Do/D a (open squares) increases strongly 
when approaching the transition. Dash-dotted line: semi em- 
prical expression by Lionberger and Russel for hard spheres. 
Solid circles: elastic spring constant k = ksT/S 2 from DWS. 
Solid line is calculated from Eq.© with ($ c = 0.625, a = 
R b /R = 0.8) 

tures 27°C and below, where the system has entered a 
viscoelastic solid phase, the swollen particles are trapped 
and the assembly is dynamically arrested. We know that 
in the swollen state the particle's mass is primarily con- 
centrated in a sphere of size Rb, whereas the dynamic 
properties are governed by the effective size R which in- 
cludes a compliant corona, R ~ Rb + L$. We can there- 
fore consider the thermal motion of the particle core, ra- 
dius Rb, as a probe of the particle interactions defined by 
the quenched particle coronas. The particle MSD shown 
in Fig. [2] clearly indicates a plateau value at long times. 
From the equipartition of energy a local spring constant 
can be defined by the maximum value S 2 of the MSD in 
the plateau region, k ~ kT/S 2 , which can be directly re- 
lated to the interaction potential ip between two spheres 
k = d 2 i/j/dr 2 . In Fig. [3] we show the temperature de- 
pendence of k in the arrested state, determined from the 
mean square displacement at t ~ 0.1s. 

As outlined above we propose that the particle corona 
can be modeled as a polymer brush of thickness L [ll|, 
12[ that extends into the good solvent from a reference 
surface located at Rb- Since typically Lq <C Rb curvature 
effects on the brush configuration and the particle inter- 
action potential should be small. Thus the force between 
two particles separated by r = 2Rb + d may be written 



in the Derjaguin approximation 24, 25] as 



F(d) = irRb / f(x)dx 



(1) 



where /(d) is the force per unit area between two flat 
surfaces 27J. The Alexander - de Gennes scaling model 
11] predicts an f(d) of the form 



f(d) ~ k B T/s 3 \(2L /d) 9 / 4 - (d/2L f/ 4 



(2) 



with an unknown prefactor of order 1. The first term 
in Eq [3] is set by the osmotic repulsion of a semidi- 
lute solution of equal density whereas the second term 
is the chain tension due to the entropic penalty upon 
stretching the grafted chains. The spring constant is 
thus k — dF(d)/dd = irRbf(d). Defining a charac- 
teristic ratio a — Rb/R < 1 we can replace d/2L 
by (r/R-2a)/{2 -2a) We can now write the 

spring constant in terms of the reduced volume fraction 
<| = $ e ^/$ c since r/R ~ 2$~ 1 / 3 . For random close 
packing conditions, $ c ~ 0.64. Note that at initial con- 
tact, r — 2R, we recover the liquid-to-solid transition. 
For simplicity we neglect the weak temperature depen- 
dence of Rb and find : 
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The spring constant directly sets the scale for the 
macroscopic shear modulus via irRG p ~ k. This can 
be obtained at a scaling level by equating the displace- 
ment energy 5 2 RG P with the thermal energy kT. For 
a multiparticle system, the same scaling relation can 
also be derived from the statistical mechanical theory 
of Zwanzig and Mountain 28|, 2t| . Figures |3] and 0] show 
the excellent agreement of Eq. [3] with the DWS-data 
using Rb/R = 0.8 and $ c = 0.625. To cover an even 
larger range of effective volume fractions we include in 
our analysis literature data of two well characterized mi- 
crogel systems The two systems are similar to ours 
except for the lower amount of BIS cross-linker used for 
one of the systems. Interestingly the moduli of the highly 
cross-linked particles and our data can be collapsed on 
single curve, whereas the particles with a lower amount 
of cross-linker are substantially softer. Moreover, both 
the more rigid and the soft particle rheology can be de- 
scribed by our model. Note that, over a limited range of 
densities, G p (& e ff) resembles a power-law with an expo- 
nent m that depends strongly on a — Rb/R. 

In this work we have focused our attention on the gen- 
eral scaling concepts that govern the linear elasticity in 
microgel pastes. We showed that a simple polymer brush 
model captures the essential physics of the interparti- 
cle interactions, regardless of the various approximations 
made. Moreover, we were able to make detailed predic- 
tions for the bulk elastic modulus over a large range of 
effective volume fractions. Our results explain the ori- 
gin of the power-law scaling of the elastic modus in these 
systems, which emerges naturally as a consequence of in- 
terparticle interactions mediated by brushlike coronas of 
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shear modulus G [Pa 




FIG. 4: Shear modulus of quenched microgels as a function 
packing density. Full symbols are derived from the diffusing 
wave spectroscopy data shown in Fig. f2]and|3] Open symbols: 
data reproduced from Ref. [j. Squares (M-5.5/115): cross- 
linker concentration 5.5 mol % BIS and hydrodynamic radius 
R H (25°C) = 115nm. Open triangles (M-l.4/141): 1.4 mol 
% BIS and R H {25°C) = 141nm. Solid lines are calculated 
from Eq.© with (<E> C = 0.625, a = R b /R = 0.8) and ($ c = 
0.65, a — Rb/R = 0.63). Dashed straight lines indicate a 
power-law scaling G v oc "3?^}/ with m = 1 + n/3. 



the microgel particle constituents. This simple picture 
should have broad utility for understanding the elastic 
properties of this important class of disordered soft solids. 
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